Possible X-ray diagnostic for jet/disk dominance in Type 1 AGN 
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ABSTRACT 

Using Rossi X-ray Timing Explorer Seyfert 1 and 1.2 data spanning 9 years, 
we study correlations between X-ray spectral features. The sample consists of 
350 time-resolved spectra from 12 Seyfert 1 and 1.2 galaxies. Each spectrum 
is fitted to a model with an intrinsic powerlaw X-ray spectrum produced close 
to the central black hole that is reprocessed and absorbed by material around 
the black hole. To test the robustness of our results, we performed Monte Carlo 
simulations of the spectral sample. We find a complex relationship between the 
iron line equivalent width (EW) and the underlying power law index (T). The 
data reveal a correlation between T and EW which turns over at T < 2, but 
finds a weak anti-correlation for steeper photon indices. We propose that this 
relationship is driven by dilution of a disk spectrum (which includes the narrow 
iron line) by a beamed jet component and, hence, could be used as a diagnostic 
of jet-dominance. In addition, our sample shows a strong correlation between 
R and T, but we find that it is likely the result of modeling degeneracies. We 
also see the X-ray Baldwin effect (an anti-correlation between the 2-10 keV X- 
ray luminosity and EW) for the sample as a whole, but not for the individual 
galaxies and galaxy types. 
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Introduction 



Time-resolved X-ray spectroscopy studies of active galactic nuclei (AGN) offer the op- 
portunity to investigate emission regions near the central black hole. In fact, X-ray spec- 
troscopy offers the clearest view of processes occurring very close to the black hole itself, 
probing matter to its final plunge into the black hole. Armed with such information, we can 
unlock the structure of the innermost regions of AGN. 

Typical X-ray spectra of AGN show an underlying powerlaw produced near the central 
black hole with signatures of reprocessed photons often present. These reprocessed photons 
show up as an Fe Ka line at ~6.4 keV and a "reflection hump" which starts to dominate near 
10 keV and is produced by the combined effects of photoelectric absorption and Compton 
downscattering in optically-thin cold matter irradiated by the hard X-ray continuum. The 
Fe Ka line has been observed in both type 1 (unabsorbed) and type 2 (absorbed) Seyfert 
galaxies. It has been attributed to either the broa d line region, the accretion disk, the molec- 
ular torus of unification models (jAntonuccilll993l ). or some combination of these. Signatures 
of reflection have also been observed in both Seyfert 1 and 2 galaxies. 

If the unification models are correct, we should see similar spectral correlations be- 
tween Seyfert 1 and 2 galaxies, with any differences easily attributable to our viewing 
angle. Regardless of the accuracy of the reflection models, we expect changes in the un- 
derlying continuum to drive changes in the reprocessing features. However, results from 
X-ray spectral studies of A GN have so far produced puzzlin g results. Samples of Seyfert 



1 obser vations from ASCA (IWeaver. Ge 



plorer (IMarkowitz. Edelson. fe Vaughan 



jord fe Yaqoobll200ll ) and Rossi X-ray Timing Ex- 



20031 ) have shown no obvious relationship between 



changes in the continuum and iron line. Several galaxies have shown an anticorrelation 
between reflection a nd /or iron line eq uivalent width and th e source flux; e.g. NGC 5548 



(Chiang et al 



2002 



2000), MCG -6-30- 15 (IPapadakis et al.ll2002h. NGC 4051 (IPapadakis et al. 



Wang et al.lll999h . NGC 5506 JPapadakis et al. 



2002 



Lamer. Uttlev fe McHardvll2000[ ). 



Recent data from Suzaku on MCG —6-30-15, on the other hand, show that the iron line and 



reflec ti on remain relatively constant while the powerlaw is highly variable (IMiniutti et al. 



20061 ). IZdziarski. Lubihski fe Smith! (119991 ) found that Seyfert galaxies and X-ray binaries 
show a correlation between the continuum slope and reflection fraction, so those with soft 
intrinsic spectra show stronger reflection than those with hard spect ra. However, othe r stud- 
ies have found e ither a shallower relationship than Zdziarski et al. (IPerola et al.ll2002l ) or an 
anticorrelation (IPapadakis et al.ll2002uLamer. Uttlev fe McHardviboOoh . 



Here we present the first results of a larger study of the X-ray spectral properties of 
Seyfert galaxies observed by the Rossi X-ray Timing Explorer (RXTE). Our full study 
consists of observations of 30 galaxies. In this letter, we focus on the spectral results from 
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the subset of 12 Seyfert 1 and 1.2 galaxies. In § [2] we present our method of data analysis, 
including our sample selection criteria (§ 12.11) . a description of our data pipeline (§ 12. 2p . and 
results of our spectral analysis (§ 12.31) . We discuss the implications of our results in § Eland 
detail our conclusions in § SI 

2. Data Analysis 
2.1. The Sample 

The RXTE public archival represents one of the largest collections of X-ray data for 
AGN, with pointed observations of over 100 AGN spanning 10 years. The RXTE bandpass 
allows the study of absorption and iron line properties of AGN spectra, as well as a glimpse 
at the Compton reflection hump. We use data from the RXTE proportional counter array 
(PCA), which is sensitive to energies from 2 to 60 keV and consists of five Proportional 
Counter Units (PCUs). Most of the sources in our sample do not show significant counts in 
the RXTE Hard Energy X-ray Timing Experiment (HEXTE), so we do not include HEXTE 
data in this study. 

To focus this study, we choose only Seyfert galaxies for which the RXTE public archive 
contained a minimum of two pointings separated by at least two weeks. We further required 
the total observed time be > 40 ks. These selection criteria led to a sample of 40 Seyfert 
galaxies. For the analysis presented here, we examine the 18 Seyfert 1 and 1.2 galaxies. Six 
galaxies were eliminated after they were put through our data pipeline (see § 12.21 for more), 
so the final sample presented here consists of 12 galaxies, listed in Table [TJ Because the data 
come from the public archive, the sample is not uniform from galaxy to galaxy or even from 
observation to observation; however, we use the Standard 2 data, which provides a standard 
data mode for these diverse observations. 



2.2. Data Pipeline 

To ensure consistent data reduction of the large volume of data, we developed a data 
pipeline. The Standard 2 data for each observation was reduced using a combination of 
FTOOLs and the Python® scripting language. The pipeline produces time-resolved spectra, 



1 Hosted by the High Energy Astrophysics Science Archive and Research Center (HE AS ARC; 
http : / /heasarc . gsf c . nasa. gov/) 
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each with a minimum of 125,000 net ph otons, which are ex tracted using standard PCA 
selection criteria and background models (IJahoda et al.ll2006l ). Sources which did not have 
sufficient net photons for even one spectrum were eliminated from the final sample (Table [T] 
shows the final sample with the 6 eliminated sources listed in the table notes). Each spectrum 
includes 1% systematic errors. We are confident in the instrument response and background 
models up to energies of ~25 keV, so we ignore channels with higher energies. 



2.3. Spectral Fitting and Results 

The data pipeline produced 350 spectra for the 12 galaxies in our sample. Each spectrum 
was fitted from 3 to 25 keV wi th an absorbed Comptqn refle ction model plus a Gaussian iron 



line. In xspec, the PEXRAV (IMagdziarz fc Zdziarskilll995l ) model simulates the effects of an 



exponentially cut-off powerlaw reflected by neutral matter and has seven model parameters: 
photon index of the intrinsic underlying power-law (r), the cutoff energy of the power law 
in keV (E c ), the relative amount of reflection (R), the redshift (z), the abundance of heavy 
elements in solar units (Z), the disk inclination angle (i), and the photon flux of the power 
law at 1 keV in the observer's frame (A). The relative amount of reflection is normalized to 
1 for the case of an isotropic source above a disk of neutral material (Q = 2ir). Adding a 
Gaussian line (energy in keV (Ep e ), physical width (er) in keV, and normalization in units of 
photons cm -2 s^ 1 ) and an absorbing column (Njj, in cm -2 ) yields a total of 12 parameters. 

We fixed the following values in PEXRAV: E c = 500 keV, Z = 1.0, and cosi = 0.95. 
This inclination represents an almost face-on disk; however, since we are seeking trends in 
the spectral parameters, rather than absolute values, the precise value is not important to 
this study. In addition, z is fixed at the appropriate value from the NASA Extragalactic 
Database for each galaxjE). After fitting all spectra to this model, we derived the mean 
Gaussian width for each source (Tabled]), then held a fixed for a second fit to the model. 
Our final model has free parameters: T, R, A, Ep e , iron line normalization and Nh- To 
prevent xspec from pursuing unphysical values of the parameters, we set the following hard 
limits: < T < 5, < R < 5, 5.5 < E Fe < 7.5 keV, and < a < 1.5 keV (for the free-cx 
fits). 

Looking at the iron line equivalent width (EW) and T, we find a complex relationship 
with a "hump" peaking near V ~ 2.0 (Figure QJi). The EW-T plot shows a correlation for 
T < 2.0 and an anti-correlation for V > 2.0, with a peak near V ~ 2.0 with EW ~ 250 
eV. We also find a strong correlation between R and V (Figured), with a best-fit line of 



^http : //nedwww. ipac . caltech . edu/ 
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R = -0.87 + 0.54 T (x 2 = 506/349 = 1.46). 

We performed a Monte Carlo simulation to determine if our results were an artifact 
of modeling degeneracies. Each spectrum in the Monte Carlo sample was simulated with 
N H =10 22 cm" 2 , r=2.0, #=1.0, E Fe =6A keV, and a Fe =0.23 keV. The flux and exposure 
times were randomly varied for each spectrum. The flux was varied by randomly choosing A 
from a uniform distribution between 0.004 and 0.06 photons keV~l cm~2 s~l. The exposure 
time was randomly generated from a uniform distribution between 300 and 11000s. The 
ranges for A and the exposure time represent the range of A and exposure for the spectra 
in the full sample. 

We generated 200 spectra: 100 simulated using an RXTE Epoch 3 response, 50 using an 
Epoch 4 response, and 50 using an Epoch 5 response, roughly corresponding to our RXTE 
sample. Each spectrum was then fitted to the same model as our full sample. The R over V 
plot (Figure [2b) clearly shows a strong correlation with a best-fit line of R = — 7.3 + 4.1 T 
(x 2 = 28.96/159 = 0.182), which strongly suggests that the observed R-T correlation is a 
result of modeling degeneracies. The correlation shows a much steeper relationship than the 
Seyfert 1 data, due to the large number of Seyfert 1 spectra showing R ~ 0. 

EW and T, however, do not suffer the same degeneracies, which is clear from the Monte 
Carlo results (Figure [lb). Based on the lack of correlation in our Monte Carlo results, we 
are confident that the shape of the EW-T plot for the data sample is real. 

To further examine the EW-T relationship, we reproduced the EW-T plot to show the 
contribution from each galaxy (Figure [3]). The radio- loud galaxies form the rising leg, with 
the quasar, 3C 273, anchoring the low T-low EW portion of the plot. The Seyfert 1 (radio 
quiet) and 1.2 galaxies tend to congregate at the peak and the falling leg of the plot. The 
one narrow-line Seyfert 1 diverges from the main cluster of points. 



X) 1 



Finally, we examined EW as a function of the intrinsic 2-10 keV X-ray luminosity (L 
using H = 70 km s _1 Mpc -1 . We fitted the data for each galaxy, each type, and the sample 
as a whole to linear and powerlaw models. The data were well-fit for either model. For 
consistency with other publications, we report here the powerlaw re sults. For the sample 



as a whole, we see an anticorrelation, i.e. the X-ray Baldwin effect (jlwasawa fc Taniguchi 
19931 ). with EW oc £-°- 14±a01 . When examining galaxy types, however, the anticorrelation 



does not always hold up (Table [T]). We find an anticorrelation in the radio loud galaxies and 
the Seyfert 1.2s, but a marginal correlation for the quasar and radio quiet Seyfert Is. 
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Discussion 



3.1. EW-Y Relationship 



The simulations of iGeorge fc Fabian! (119911 ) for the observed spectrum from an X-ray 
source illuminating a half-slab showed that the spectra should include a "Compton hump" 
and an iron line. They found that the iron line EW should decrease as the spectrum softens. 
This is easy to understand, since as the spectrum softens (Y increases), there are fewer 
photons with energies above the iron photoionization threshold. Our results show that the 
relationship between EW and T is not quite so simple. We find a correlation between EW 
and T when T < 2 and an antico rrelation when Y > 2. Other researchers hay e found a 
correlation for Seyfert 1 samples (IPerola et al.l 120021 ; iLubihski fc Zdziarskil 20011 ) , but the 



galaxies in their samples primarily fell in the T < 2 region. iPage et al.l (120041 ) also find that 
their data suggest a slight correlation for a sample of radio loud and radio quiet Type 1 
AGN. 

A close examination of our EW-Y plot shows that the data for different galaxy types 
progresses across the plot. The plot is anchored at the low-r, low- EW end by the quasar, 
3C 273, in our sample. The rising arm of the plot, T ~ 1.5 — 2.0 and EW ~ — 300 eV, is 
primarily formed by radio loud Seyfert 1 galaxies. The radio-quiet Seyfert 1 galaxies cluster 
near the Y ~ 2.0, EW ~ 300 eV peak of the hump, and the radio-quiet Seyfert 1.2 galaxies 
form the falling arm of the plot for Y > 2.0. 

Physically, the most obvious difference between these sources is the presence or absence 
of a strong jet. We propose that this relationship is driven by the degree of jet-dominance 
of the source. The iron line features are associated with the X-ray emission from the disk. 
Since the disk is essentially isotropic, it will excite an observable iron line from matter out of 
our line-of-sight. On the other hand, the jet is beamed away from the obvious configurations 
of matter in the system and, more importantly, is beamed toward us in the quasar and radio- 
loud sources. Both of these jet-related phenomena reduce the observed equivalent width of 
any iron line emission associated with the jet continuum. 

In order for the Y to increase as the jet-dominance decreases, the jet in these sources must 
have a hard X-ray component, which implies that the radio-loud Seyferts in our sample are to 
be associated with low-peaked BL La c objects (LBLs). BL Lac o bjects show two broad peaks 
in their spectral energy density plots (jGiommi fc Padovanill 19941 ). with the lower-energy peak 
due to synchrotron emission and the higher-energy peak due to inverse Compton emission. 
BL Lacs are divided into two classes, depending on where the peaks occur: high-peaked BL 
Lacs (HBLs) and LBLs. The X-ray continuum in the HBLs is rather soft, since we are seeing 
the synchrotron spectrum cutting off in these sources. LBLs, on the other hand, tend to 
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have a harder X-ray continua, since we are observ ing: well into the inverse Compton part of 
the spectrum (IDonato. Sambruna. fc Gliozzil 120051 ) . 

We also note that much of the falling arm of the EW-Y relationship is formed by MCG 
—6-30-15. Recent observations of MCG —6-30-15 by Suzaku ha ve shown that the refl ection 
component, including the iron line, remains relatively constant iMiniutti et al.l (120061 ). We 
would expect, then, that as Y increases, the EW should decrease, which is exactly what we 
see in our data. 



3.2. R-Y Relationship 

Significant degeneracies between the photon index, absorbing column, and reflection 
fraction can easily lead to false conclusions about spectral correlations. These degeneracies 
occur as these three parameters trade off against each other in the modeling process, an effect 
that is especially strong in the RXTE bandpass. Our R-Y plot shows a strong correlation 
which is mimicked in our Monte Carlo results. The few points that lie under the main 
concentration are likely to be outliers, and not indicative of a subclass of galaxy. These 
points all come from spectra that have been fitted to have Nh = 0, and are primarily 
radio-loud galaxies. We conclude that the observed R-Y correlation in our sample cannot be 
trusted as a real correlation. 



3.3. EW-L X Relationship 



Looking at the EW-L X relationship, we do see the X-ray Baldwin effect for our sample 
as a whole, with a slighly shallower an t icorre latio n than reported elsewhere. W e find EW oc 
L~ 014 , whe reas llwasawa fc Taniguchil (119931 ) and IJiang. Wang fc Wang] (2006) find EW oc 



-0.20 



and IPage et al.l (I2004T ) find EW oc L 



-0.17 



However, when IJiang. Wang Wang 



(120061 ) exclude the radio loud galaxies from their sample, they find EW oc L, 



-o.io 



We find, though, that when we examine our data on a galaxy-by-galaxy or type-by-type 
basis, the effect is not consistent from source to source. At this point, we cannot determine 
if these variations are real or are simply due to the small number of spectra for some of our 
galaxies and types. 
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4. Conclusions 

We have examined time- resolved spectra of 12 Seyfert 1 and 1.2 galaxies observed by 
RXTE over seven years. We find a complex relationship between the iron line equivalent 
width and the continuum slope, with a correlation for T < 2 that turns over to an anticor- 
relation for T > 2. We propose that this relationship is a possible diagnostic for jet- versus 
disk-dominated sources, where jet-dominated sources show a correlation between EW and 
T, and disk-dominated sources show an anticorrelation. We also see a strong correlation 
between V and R which is likely an artifact of modeling degeneracies caused by the interplay 
of T, R, and n H in the RXTE bandpass. Finally, we observe the X-ray Baldwin effect for 
the sample as a whole, but not for each galaxy and galaxy type individually. 

This research has made use of data obtained from the High Energy Astrophysics Science 
Archive Research Center (HEASARC), provided by NASA's Goddard Space Flight Center. 

This research has also made use of the NASA/IPAC Extragalactic Database (NED) 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and Space Administration. 
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Table 1. Sample of RXTE-ohserved Seyfert 1 and 1.2 galaxies 8. 
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3 
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a The following sources were eliminated after running the data pipeline described 
in the text, due to having no spectra with at least 125,000 net photons: Mkn 110, 
PG 0804+761, PG 1211+143, Mkn 79, Mkn 335, and PG 0052+251. 

b Seyfert type based on the NASA Extragalactic Database 

c Total number of spectra extracted using our data pipeline (§ 12.21) . 

d The average physical width of the Fe Ka line for all spectra from a source when 
fitted to the absorbed powerlaw model with Compton reflection and Gaussian iron 
line (§E3D. 
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e Results of fitting the X-ray luminosity over EW plot to a powerlaw model; e.g. 
EW oc L", where L x is the 2-10 keV X-ray luminosity in ergs s _1 and EW is the 
iron line equivalent width in eV. 

f One 3C 120 spectrum shows a flare, where L x jumps by ~ 6x. The number 
quoted above excludes this point from the sample. If we include the flare, we find 
EWocL° x 07i+0 - 18/ -°- 25) . 
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Fig. 1. — Iron line equivalent width in eV (EW) versus powerlaw photon index (T) for the 
Seyfert 1/1.2 sample (a) and for the Monte Carlo simulations (b). 
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Fig. 2. — Reflection fraction (R) versus powerlaw photon index (r) for the Seyfert 1/1.2 
sample (a) and for the Monte Carlo simulations (b). In both plots, the line shows the 
best-fit linear model for the Monte Carlo simulations. 
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Fig. 3. — The iron line equivalent width in eV versus the powerlaw photon index. This plot 
is similar to the left panel in Figure [U but with each galaxy plotted with a separate symbol. 
The open circles are 3C 111, open squares are 3C120, pluses (+) are 3C273, open triangles 
are 3C 382, open diamonds 3C 390.3, open stars Akn 120, open crosses Fairall 9, filled circles 
IC 4329A, filled squares MCG -6-30-15, filled triangles Mkn 509, filled stars NGC 7469, and 
asterisks (*) TON S180. 



